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Abstract 

The possibility of the decomposition of the three dimensional (3D) Gross-Pitaevskii equation 
(GPE) into a pair of coupled Schrodinger-type equations, is investigated. It is shown that, under 
suitable mathematical conditions, solutions of the 3D controlled GPE can be constructed from the 
solutions of a 2D linear Schrodinger equation (transverse component of the GPE) coupled with a 
ID nonlinear Schrodinger equation (longitudinal component of the GPE). Such a decomposition, 
called the 'controlling potential method' (CPM), allows one to cast the above solutions in the form 
of the product of the solutions of the transverse and the longitudinal components of the GPE. The 
coupling between these two equations is the functional of both the transverse and the longitudinal 
profiles. The analysis shows that the CPM is based on the variational principle that sets up a 
condition on the controlling potential well, and whose physical interpretation is given in terms of 
the minimization of the (energy) effects introduced by the control operation. 
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I. INTRODUCTION 



Since its experimental observation the study of the three dimensional (3D) collective 



and nonlinear dynamics of the Bose Einstein condensate (BEC) [2] in an external potential 
trap |3j has received a great deal of attention by a very wide scientific community and in the 
investigations concerning fundamental physics, by mathematical physics and sophisticated 
technological applications [4]. Although rapid scientific and technological advances have 
been achieved in this area, finding the exact analytical 3D solutions of the Gross-Pitaevskii 
equation (GPE) [3], that correspond to the coherent state of a BEC in a suitable external 
potential well (such as soliton-like structures), still remains a challenging task for physicists 
and mathematicians. 

n n 

A number of valuable approximative analytical [5] and numerical evaluations [6J have been 
presented in the literature and have been adequately compared with a very wide spectrum 
of experimental observations. The experience gained from these investigations may suggest 
;he idea that a BEC's dynamics exhibits the features of a nonlinear non-autonomous system 
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for which it seems to be necessary to include some control operations in order to allow 
the existence of coherent structures. In particular, to retain the 3D coherent stationary 
structures of the BEC for a long time, suitable " ad hoc" time-dependent external potentials 
and control operations are known to be necessary Furthermore, in the presence of an 
inhomogeneous time-dependent external potential one encounters some difficulties to find 
exact soliton solutions in one or more dimensions, although several kind of solitons have 
been found in certain approximations [9] . Consequently, one easily arrives to the conclusion 
that, in order to get exact soliton structures, some sort of the 'control of the system' seems 
to be necessary. This implies that the correct analysis of the system should include a control 
potential term in the GPE which is to be determined dynamically by the system itself. In 
principle, this procedure may be extended to an arbitrary 'controlled solution' with the 
appropriate choice of the external potential (so-called 'controlling potential' [h]] ). In fact, 
a controlling potential method (CPM) has been recently proposed in the literature and 
used to find multi-dimensional controlled localized solutions of the GPE. In the preliminary 
investigations ll||, this method has established reasonable experimental control operations 



that ensure the stability of the solution against relatively small errors in the experimental 
realization of the prescribed controlling potential. The main goal of the CPM is to fix the 
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type of the desired controlled solution and to find the appropriate family of the controlling 
potentials. Then, the set of suitable mathematical conditions has to be found allowing us 
to select the desired solution, with the adopted controlling potential. 

In this manuscript, we develop an analytical procedure to construct exact three dimen- 
sional solutions of a controlled Gross-Pitaevskii equation, by improving the CPM. To this 
end, we develop the theory of the BEC control based on two decomposition theorems leading 
to suitable physical conditions to express the BEC wave function as the product of a 2D wave 
function and a ID wave function, taking into account the 'transverse' and 'longitudinal' BEC 
profiles, respectively. Such a factorization allows us to decompose the 3D controlled GPE 
into a set of coupled equations, comprising a 2D linear Schrodinger equation (governing the 
evolution of the 'transverse' wave function), a ID nonlinear Schrodinger equation (governing 
the evolution of the 'longitudinal' wave function) and a variational condition involving the 
controlling potential. The requirement for the minimization of the effects introduced by the 
control operations (i. e. the requirement that the average of the controlling potential over 
the transverse plane is equal to zero) allows us to determine explicitly the self-consistent 
controlling potential which also plays the role of the coupling term between transverse and 
longitudinal BEC dynamics. 



II. CONTROLLED GROSS-PITAEVSKII EQUATION 

It is well known that the spatio-temporal evolution of the ultracold system of identical 
atoms forming a BEC in the presence of the externa! potentia! U M within the mean 
field approximation, is governed by the three dimensional Gross-Pitaevskii equation |3j, viz., 

ih— = -—V 2 ^ + NQ\y\ 2 y + U ext (r,t) (1) 

where \I/(r, t) is the wavefunction describing the BEC state, m a is the atom mass and Q is 
a coupling coefficient related to the short range scattering (s-wave) length a representing 
the interactions between atomic particles, viz., Q = 4irh 2 a/m a , and N is the number of 
atoms. Note that the short range scattering length can be either positive or negative. We 
assume that U ext is the sum of the 3D trapping potential well, U trap , that is used to confine 
the particles of a BEC, and the controlling potential U con t r which will be determined self- 
consistently. We conveniently introduce the variable s = ct (c being the speed of light) and 
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divide both sides of Eq.flT]) by m a c 2 , and we use the notation 
U e xt{T,t) _ Ut rap {r,t) U cont r(r,t) 



m a & rn a c z 
Eq. (CQ) can be cast in the form 



+ 



m a c z 



= Vtrapif, s) + V con tr(j, s) 



(2) 



(3) 



ds 2 

where x/j(r,s) = \^(r,t = s/c), Ac = h/m a c 2 is the Compton wavelength of the single atom 
of BEC and g = NQ/mc 2 . 

In this paper, we will investigate the properties of Eq. (j3J) and V^ on4r that enable the 
existence of the controlled 3D solutions in the factorized form 

^(r, s) = ip±(r±,s) V 2 (r_L,z, s) , (4) 

provided that Vt rap can be split into two parts, as 

V trap {r,s) = V ± (t ± ,s) + V z {z,s) (5) 

where, in Cartesian coordinates, r = (x, y, z) and r_|_ = (x, y) denotes, by definition, the 
'transverse' part of the particle's vector position r. We also refer to z as to the 'longitudinal' 
coordinate. 



By substituting Eqs. (Tj0) and (jSJ) in Eq. (jBJ), we easily get: 



Ac 



d^ z \ 2 c d 2 ^ z 

1 Ac n ~r 



ds 2 



i Ac + ^ViV± - V ± (r ± ,sW ± 
as 2 



\V Z + Vcontr + 9 IV'xl 2 IV'zl 2 ) i>. 
0. 



(6) 



where, in Cartesian coordinates, Vj_ = xd/dx + yd/dy. 

Let us define as 'controlled parameter' the following time-dependent quantity: 

Qid{s) = q J d 2 f ± \ip ± f ; 

and the following linear and nonlinear operators, respectively: 

\l d 2 



-yVi + ^(r ± , S ) 



2 dz< 



+ V z (z, s) + q 1D (s) \^ z (r ± , z, s)\ 2 + V 



(7) 

(8) 
(9) 



where Vq is an arbitrary real constant. Then, Eq. ([6]) can be rewritten as: 



(. v. j i Ac - H± ) < '_ ' 



i\ c -^-- H z ) 1p z + (ftjj(s) - <#±| 2 ) IV 1 * IV* + (^0 - Kontr) ^2 



X 2 



v 2 ^ + 2 • v ^ 



0. 



(10) 
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III. THE DECOMPOSITION PROPERTIES OF THE CONTROLLED GROSS- 
PITAEVSKII EQUATION 

By the definition of the controlling potential, V contr depends both on ip± and ip z . In 
particular, we assume here that the space and time dependence of V contr is given also through 
P±(r±,s) = |^±(r±,s)| 2 , viz., 



Vcontr — V con tr (p±(T±, S), Z, s) . 



(11) 



Moreover, defining also the following functional of p±: 




(12) 



the following theorem holds: 



DECOMPOSITION THEOREM 1. 



iI> z (t±,z, s) = ip z (z, s), 



(13) 



and ip±(r±, s) is the solution of the following 2D linear Schrddinger equation 




(14) 



and V is a stationary functional (with respect to variations 5p± of p±), assuming the value 
V — Vq, conditioned by the constraints 




(15) 



(normalization condition for ip±), and 




(16) 



then ip z is the solution of the following ID nonlinear Schrddinger equation 




(17) 



and V c 



contr is given by 




(18) 
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To prove this theorem, first of all, we note that the assumptions (1131) and (fT4j) allow us to 
reduce Eq. (TTU1) to 



(i Ac ^ - + - <#J 2 ] |V»*|V* + (^0 - Vcontr) fa = . (19) 

Secondly, the required stationarity of V with respect to variations 5pj_ of p_i_ implies that 

5V + a(z,s)5 J p±d 2 r ± + p(z,s)6 J pld 2 r L = 0, (20) 

where a(z, s) and f3(z, s) are Lagrangian multipliers. Taking into account Eq. ffl2l . condition 
(|20l allows us to solve the corresponding ordinary inhomogeneous first-order differential 
equation for V contr where p± plays the role of the independent variable and z and s are 
parameters, yielding the following general solution 

hi z s) 

V cont r(r±, z, s) = — - J — r - a(z, s) - /3(z, s)p_ L (r_ Ll s) , (21) 
P±(r±,s) 

where h(z, s) is an arbitrary function. Actually, to ensure the convergence of the integral in 
the definition of the functional V, see Eq. (Tl2l) . it is easy to see that we must have h(z, s) = 0. 
Consequently, the appropriate V contr satisfying the stationarity condition V = Vq is given by 



[r ± ,z, s) 



Q 

which after the substitution in Eq. ffl9l gives 



P(z,s) + V , (22) 



(i A c ^ - H^J ij z + [q 1D (s) - q\4> ± \ 2 ] (|^| 2 - P/q) fa = . 



(23) 



Now, according to the hypothesis f|T3|) . to preserve the r^-independence of ip z , Eq. (j2"3"j) can 
be satisfied only when 

P(z,s)=q\fa(z,s)\ 2 , (24) 
which immediately implies that Eqs. ( fl7|) and ( fl8|) are satisfied. 

DECOMPOSITION THEOREM 2. 

Lei its suppose that ip z = tp z (z,s) is the solution of the ID nonlinear Schrddinger equation 
Then, the functional V given by [T2}) and conditioned by the constraints ( T73j) . and (T76j) ; 
zs stationary (with respect to variations 5p± of p±), V — Vq if, and only if, t/j± = ip±(r±,s) 
is the solution of the 2D linear Schrddinger equation [Lfy . 
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To prove this proposition, we observe that since ip z (z, s) satisfies Eq. (iT7|) . Eq. (fTUj) becomes 
(i A c ^ - #±J V.L + [(?u>00 - <#±| 2 ) \^z\ 2 + (V - V contr )] ij> ± = 0. (25) 

By multiplying the latter on the left by ip± and integrating over all the transverse plane, we 
easily obtain 

/ ^( 2A ^-#±)^V + ^-V[p ± ;z,s]=0, ( 26 ) 

where constraints ( fT5l) and ( TlBl) have been used. Consequently, if ip± satisfies Eq. ( 1X4*1) . then 
V is a stationary functional with the value V = Vo, conditioned by (|15|) and ffTB]) . Conversely, 
the assumed stationarity of V implies that the functional form of V contr with respect to r_i_, 
z and s is given by Eq. ([22"]) . which substituted in Eq. (|25|) gives 

f i Ac ^ - H ± j ij, ± + (q 1D (s) - q\^\ 2 ) (|V,| 2 - /?(*, = . (27) 

However, if /3(z,s)/q ^ \ip z (z, s)\ 2 , then ip z would be also function of rj_ which would con- 
tradict the assumption tp z = ip z (z,s). It follows that j3(z,s)/q = \ip z (z, s)\ 2 and, in turn, 
that Eq. ( TH|) is satisfied. 

The results presented above allow us to draw the following conclusion. 

If "0_i_( r _L, s) and ip z (z, s), are two complex functions which are exact solutions of the 2D linear 
Schrddinger equation (fl^p and the ID nonlinear Schrddinger equation [TJ\ ), respectively, 
provided that V cont r is given by Eq. ( TJg|) . the function ip(r,s) = ij;±(r±,s)ip z (z,s) is the 
exact solution of the controlled 3D Gross-Pitaevskii equation |3)]. 

Of course, the inverse is not necessarily true. In fact, it is easy to see that, in principle, 
it is not true that an arbitrary solution of Eq. (J3j) can be expressed as the product of two 
wave functions ^±( r -L) s ) an d ?p z (z,s) that obey the Eqs. ( 1141) and ( |T7I) . respectively. In 
other words, we can decompose the controlled 3D GPE ([3D into the system of equations 
(|T4|) . (|T7|) and (|T8j) only for the subset of its solutions of the type (jlj). However, using such a 
decomposition we are able to solve Eq. ([3]) and to obtain a wide spectrum of exact solutions 
of the type (j4]). 



IV. CONCLUSIONS AND REMARKS 



In this paper, we have presented some mathematical properties of the controlled 3D 
GPE ([3]). After formulating and proving two decomposition theorems, we have found the 
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mathematical conditions that make possible the construction of the solution in a factor- 
ized form, i.e. ip(r,s) = ip±(r±,s) ip z (z,s), where ip±(r±,s) and ip z (z,s) satisfy the 2D 
linear Schrodinger equation (i \ c d'ip±/ds = H±i{j±) and the nonlinear controlled nonlinear 
Schrodinger equation (i \ c dip z /ds = H z ip z ), respectively. The results presented here improve 
the formulation of the recently proposed Controlling Potential Method [lol. 

It is worthy observing that the set of equations (fill) . (TTTT) and (|T8|) opens up the possibility 
to find the controlled solutions of the type (jlj) which exhibit the quantum character in 
the transverse part (superposition principle with consequent interference effects) and the 
classical character in the longitudinal part (due to the nonlinearity of the ID nonlinear 
Schrodinger equation), although the entire solution of the controlled 3D GPE is nonlinear 
and, therefore, has a classical character. By means of suitable controlling and trapping 
potentials, this possibility would allow, for instance, for a very stable soliton-like longitudinal 
profile of the BEC whose transverse profile would have a quantum character as a result of 
the quantum interference at the macroscopic level. 

Note that, when ip± satisfies Eq. (tHj), according to definition (T121) . V represents the 
average of V con t r in the transverse plane. The value of this average corresponds to the 
arbitrary constant Vq. Without loss of generality, we put Vq = 0, viz. 



This way, among all possible choices of V contr , we adopt the one which does not change the 
mean energy of the system (note that the average of the Hamiltonian operator in Eq. (J3]) is 
the same with or without V contr ) and thus minimizes the effects introduced by our control 
operation. 

In our forthcoming papers, we will use the method developed in the present paper to solve 
exactly the 3D controlled GPE with a 3D parabolic potential trap. We find the controlled 
envelope solutions in the form of localized as well as periodic structures for which suitable 
stability analysis is performed. 
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